Repulsive Casimir Force in Chiral Metamaterials 
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We demonstrate theoretically that one can obtain repulsive Casimir forces and stable nanolevitations by using 
chiral metamaterials. By extending the Lifshitz theory to treat chiral metamaterials, we find that a repulsive force 
and a minimum of the interaction energy exist for strong chirality, under realistic frequency dependencies and 
correct limiting values (for zero and infinite frequencies) of the permittivity, permeability, and chiral coefficients. 
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Following the original Casimir paper for the attraction 
of two media, 1 and 2 occupying half spaces, z < and z > d, 
respectively, and such that the electromagnetic fields are con- 
fined exclusively in the vacuum region between them, Lifshitz 
generalized the calculation of this force to the case that 
these two media are characterized by frequency-dependent di- 
electric functions ei(ui) and £2(w). Subsequently, there was 
further generalization to general bi-anisotropic media 01 . The 
formula for the force or the interaction energy per unit area 
can be expressed in terms of the reflection amplitudes, r°; b 
(j = 1, 2) 01, at the vacuum/medium j interface, giving the 
ratio of the reflected EM wave of polarization a by the in- 
coming wave of polarization b. Each a and b stands for either 
electric (TM or p) or magnetic (TE or s) waves. The frequency 
integration is completed along the imaginary axis by setting 
co = i£. The formula for the interaction energy per unit area 
becomes 
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For isotropic media, the off-diagonal terms in Eq. (O vanish 
and 
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where Kj = wky + e^^/c 2 and fij is the permeability of 



medium j. 

In most cases the resulting Casimir force between the two 
media separated by a vacuum region is attractive. There is in- 
creased interest recently |@, Q, & 0] in determining whether 
there is a combination of media 1 and 2 capable of producing 
a repulsive force. There have been mainly three mechanisms 
to obtain repulsion for the Casimir force: (1) Dzyaloshinskii's 
Casimir repulsion 10]: Immersing the interacting plates of ei 



and 62 in a fluid of 63 and, moreover, satisfying the condi- 
tion ei(i£) < £3(«£) < £2(i£); (2) Boyer's Casimir repulsion 
|01 : Based on an asymmetric setup of mainly (purely) non- 
magnetic/vacuum/mainly (purely) magnetic; (3) Leonhardt's 
Casimir repulsion [g]: Employing a perfect lens sandwiched 
between the interacting plates. The possibility for a transi- 
tion from an attractive to a repulsive force as the distance 
d decreases (corresponding to a minimum of the interaction 
energy) leads to nanolevitations and opens up many oppor- 
tunities for application, e.g., almost frictionless operation of 
nanomotors. Even through Capasso's group experimentally 
realized the repulsion, based on the theoretical prediction of 
Dzyaloshinskii et al. |6], this kind of system still has friction 
because of the existence of the liquid. Leonhardt's Casimir 
repulsion needs a perfect lens with simultaneously negative 
dielectric permittivity and magnetic permeability, which are 
extremely difficult to obtain at optical wavelengths. Finally, 
Boyer's Casimir repulsion proposal faces the essential ob- 
stacle that such nontrivial magnetic materials in the optical 
regime do not exist in nature, and, therefore, it relies on the 
nontrivial possibility of developing new artificial negative in- 
dex metamaterials (NIMs). 

In this letter, we examined realistic non-chiral metamate- 
rials and we concluded they do not give a repulsive Casimir 
force. However, we found that chiral metamaterials are excel- 
lent candidates to realize the repulsive Casimir force. The ex- 
istence of a repulsive Casimir force depends upon the strength 
of the chirality. We present analytical arguments that strong 
chirality gives a repulsive force, supported by numerical cal- 
culations. 

Negative index metamaterials ifioll . because of their reso- 
nance magnetic response, offer more flexibility and, hence, 
more promise for achieving a repulsive Casimir force, based 
on Boyer's prediction. Indeed, in recent papers, Rosa et al. 
1 1 ill found a repulsive force in a range of values of d for a 
combination of a mainly nonmagnetic Drude-modeled silver 
and a magnetic NIM. This result was obtained 1 1 ill through 
the employment of a Lorentz type of magnetic permeability 
of the form = 1 — O 2 / (lo 2 — uj^ + ijoj). This form pro- 
vides the opportunity to use an f2 large enough as to satisfy 
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the condition > e(i£) and obtain thus Boyer's Casimir 

repulsion. For the reasons stated below we consider a Lorentz 
type frequency dependence of fi(uj) unphysical. Instead we 
employed the following realistic expression for n{u>): 
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where \a\ is usually much smaller than one and A = a in 
order to satisfy the physical requirement that (i(u>) — > 1 as 
lo — ► oo. It must be stressed that the realistic expression ©, 
although almost identical to the Lorentz form for £j around the 
resonant value uj m , produces radically different results than 
the Lorentz one as far as the Casimir attraction is concerned. 
As the authors of Ref. have found out (and we have con- 
firmed), expression (|4]i (with a = and A / 0) combined 
with the form of Eq. © below for e(u>) does not produce 
repulsion. This is also true for the realistic case of a = A. 

The uj 2 dependence of the numerator of the resonance term 
follows from the equivalent circuit approach 11211 and from the 
Maxwell's equations in the low frequency regime as stated in 
Ref. fllll . It is confirmed by the retrieval procedure in actual 
SRR based and fishnet metamaterials. Of course, it is possi- 
ble to have more than one resonance term in Eq. but their 
coefficients must satisfy the relation a — J2i At = to obtain 
the correct limiting value of /i(oo) = 1. Besides cases having 
a = A, we also examine the case a — and A ^ (which 
produces the incorrect limiting behavior, fi(u>) = 1 — A as 
u> — > oo). The reason for this unphysical choice is to deter- 
mine the role of the uj = oo value of fx(u>). The most general 
form of the frequency-dependence of the dielectric function is 
the sum of the Drude term and several Lorentz-type resonance 
terms. If only one resonance term is kept, we have 
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We have calculated the Casimir force using for material 1 
and material 2, ei, \i\ and 62, M2, as in Eqs. (|4]i and (0 with 



u) m , OJp 
,2 



(including 



0, u:i ± 



several values of A, 

0, and ojp[ ^ 0, uj 2 = 0). Among these values, we in- 
cluded realistic values as they were obtained by our retrieval 
approach in various fabricated and/or simulated NIMs. The 
Casimir force turned out to be attractive in all cases we calcu- 
lated. See the triangle and diamond curves in Fig. Q] 

Recently, a lot of experimental work on chiral metamateri- 
als (CMMs) fabricated by planar technologies have been pub- 
lished JI3I1 . For such artificial materials, the constitutive equa- 
tions have the form 
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where the coefficient n has the following frequency depen- 
dence for the chiral metamaterials 11411 : 
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which is the same as Condon model for homogeneous chiral 
molecular media Il5ll . 

For such CMMs, the reflection elements can be expressed 
as follows, assuming the electromagnetic wave is from vac- 
uum to chiral metamaterials Big], 
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where, K± 



v /k[T^|iv^, n± (it) = V^MhhW) ± 

and (J.j(i£) are the relative permittivity and permeability of 
the plate j, respectively, and Kj is the chirality coefficient. 
Although n± are complex, the reflection elements, r's, are still 
purely real because \+ = X- ■ 
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FIG. 1: (Color online) Casimir interaction energy per unit area E/A 
(in units of hckg) versus kod; ko = lor/c. The triangle curve 
corresponds to a — A — 0.001, k — (no chirality), uj m = 
ujri 7m = 7.R = 0.05wfl, Upi =0, Lo e — lur for material 1, while 
a = A — 0, u)pi = Wlur, j p i — 0.05uj p i, oj e = for material 2. 
The diamond curve is the case with a = A = 0.001, k = 0, ui m = 
lor, = jr = 0.05lor, u)pi = 0, w e = ujr. The squares curve 
is the case with a — A = 0.001, to K i — lo K 2 = Q.Qujr, u) m = 
lo k r = ojr, 7 m = -y K = 7a = 0.05ujr, ujpi — 0, cu e = lor- Fi- 
nally, the circle curve shows repulsion for kod < 0.0586 and a stable 
equilibrium point at kod = 0.0586; the parameters are the same as 
for the square curve except for co K i — lo K 2 = 0.7lor. 

Here, we consider first a special setup with two identi- 
cal chiral metamaterial plates with the following parameters: 
ei = £2 = e; Mi = M2 = MS K i = K 2 = K. We suspect that the 
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chirality coefficient, k, may provide sufficient new freedom to 
drive the force to negative values (i.e., repulsive) at least for 
some range of values of d. From Eq. ([1) it follows Ref. 1 1 ill 
that a negative value of the Casimir force is favored by mak- 
ing the quantity / = Tr[D _1 (l — D)] as negative as possible 
over as broad a range as possible of the parameters and the 
integration variables. This quantity, 7, has the same sign as the 
quantity F given below: 
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Because r sp is purely real as shown in Eq. (0, it is clear 
from Eq. (0 that the chirality by introducing the off-diagonal 
quantity r sp provides the possibility, for large enough r sp , to 
make the numerator in Eq. © negative, while keeping the 
denominator positive. Thus, the chirality, if strong enough, 
is expected to lead to a repulsive Casimir force. This expec- 
tation is confirmed by the numerical evaluation of the inter- 
action energy per unit area as shown in Fig. Q] Indeed, for 
large enough chirality parameter, uj k \ = lo k2 = 0.7ojr, we 
have a very interesting situation of an attractive force in the 
range d > do (where in the present case do = 0.0586c/wr) 
and a repulsive case for d < do. Thus, a stable equilibrium 
distance emerges, d = do, reminiscent of the bond length in a 
diatomic molecule. There is a critical value of u K , io K = lo c k , 
such that for lo k < lo c k there is no repulsive regime for any 
value of d, while for lo k > u) c K , there is a distance do, a func- 
tion of uj K , do(ui K ), such for d < do(io K ) the force is repulsive. 
For the numerical values used in our present case, the critical 
value of lo c k is equal to w° K = 0.612u;# for a(= A) = 0. As 
shown in Fig. |2(a)| the critical value io c K is a function of the a 
with its minimum value io c K = 0.607 obtained for a ~ —0.09. 
Furthermore, the relation do versus lo k (for uj k > u>%) is an 
increasing almost linear function of u) K , as shown in Fig. |2(b)| 
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FIG. 2: (Color online) (a) The critical value of chirality lo% versus 
a(— A), for two identical CMM plates, (b) The equilibrium distance 
fcodo versus lo k for a = A = 1CP 3 . For lo k > io% = 0.612lj_r, the 
value of the equilibrium distance kodo corresponds to the minimum 
of the energy as shown by the open circle curve in Fig. [TJ 

The question raised by the present novel approach to a pos- 
sible repulsive Casimir force is whether real chiral metamate- 
rials can be fabricated with ui K larger than the critical one lo c k . 
Our own chiral metamaterial presented in Ref. H 1 3fl has an 
co K ~ 0.3w£. However, this metamaterial was designed and 
fabricated before the critical importance of chirality for sta- 



ble Casimir nanolevitation was even suspected; thus, there is 
room for new designs to raise the value of io K possibly above 
the critical value. We are currently working on this theme. We 
don't know whether or not general physical considerations re- 
strict the size of the chirality factor ui K and thus we cannot 
be sure whether the critical value of lo l k is reachable. Models 
based on a single l oop (see the books of Lindell et al. Hill and 
Serdyukov et al. 11711 ) produce a relation between the elec- 
tric, a ee , the magnetic, a mm , and the cross polarizabilities, 
a em , a me : a ee a mm = a em a me . This relation, valid when 
u m = u R = uj kR , 7 m = 7 K = 7 K and uj pl = a = 0, A^0, 
shows that the critical value uj c k is almost reachable under the 
optimum condition A — ► u>^/u> R . 
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FIG. 3: (Color online) Casimir interaction energy per unit area E/A 
(in units of hck^) versus hod of the two identical CMM plates config- 
uration for different chiral strengths oj k 's. a = 0, A = 0.2, lu p i = 

0, UJ m = LO K R = LOR, Je = 7m = Jk = O.Q5ujR, LO e = UR- 

In Fig. [3]we present results for the energy per unit area E/A 
versus the dimensionless distance kod for chiral metamateri- 
als with n(ui) given by Eq. © with a = and A = 0.2. 
We repeat here this choice violates the physical requirement 
of fJb{uj) — * 1 as to — > oo. Nevertheless, we present these 
results here in order to show that unphysical frequency de- 
pendence of the response functions may produce the resulting 
behavior which is qualitatively different from that presented 
in Fig. Q]in the sense that now two equilibrium points, d\ and 
cfe (d\ < d 2 ), may appear, the first is unstable equilibrium 
and the second is stable equilibrium. Furthermore, one cannot 
exclude the possibility that a more complicated fJ,(ui) satisfy- 
ing the condition yu(oo) = 1 and producing results as those 
in Fig. [3] may exist. In spite of this unphysical behavior of 
jti(oo) (ju(oo) = 0.8 instead of jtt(oo) = 1), one expects to 
produce no repulsive force if regular metamaterials (with no 
chirality) are employed. The reason is that both of the inter- 
acting plates are mainly nonmagnetic with /Lt(i£) < 1 < e(i£) 
at all frequencies. For a not so large chirality (lj k = 0.94 
(circles)), one can easily see from Fig. [3] there is only an 
attractive Casimir force for all distances. However, as chi- 
rality increases, (u> K — 0.96 (squares)), the energy tries to 
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develop a minimum, and the Casimir force corresponding to 
the slope of the energy becomes smaller. At u K = 0.98 (dia- 
monds), there is a minimum at k d 2 — 0.21 and a maximum 
at k di ~ 0.09; the peak value of the energy is less than zero, 
the energy value at k$d = oo. If chirality increases further, 
oj k = 1.00 (triangles), the sign of the energy is reversed and 
becomes positive in a certain range. This is an interesting case 
that gives a repulsive Casimir force within a range of distances 
between d\ and d2 . It forms a potential barrier to block the 
two interacting plates sticking to each other. Similar results to 
those in Fig. [3]were also obtained for the case where we used 
e(u>) = 2 — uj 2 /(lo 2 — uj 2 + This frequency-dependence 
is obtained experimentally ft 1 Oil for realistic metamaterials, but 
only close to the resonance behavior; such a dependence ex- 
tended to u> — > oo violates the condition of e(uj) — > 1 as 

LU — > OO. 

In discussing these results we must keep in mind that for 
kod <C 1 the main contribution to the integral in Eq. (fl]i comes 
from large £ and fcii values with the ratio feii/£ 3> 1, as argued 
by Landau et al. [18] and confirmed by our numerical calcu- 
lations. Under these conditions k ~ fcii and the integrand in 
Eq. (l} takes the form /(£, e~ 2fc n d ). By setting x = 2k\\d, 
it follows immediately from Eq. ([TJ that E(d)/A oc d~ 2 
and F(d)/A oc d~ 3 ; the contribution of the chiral term to 
/(£, e~ x ) is negative and, thus, for large enough chirality the 
force in the d — > limit becomes repulsive. On the other 
hand, in the opposite limit d — > oo, because of the factor 
e~ 2Kd , the main contribution to the integral comes from the 
range < £ < (c/d) and < fc|| < (d^ 1 ), where the inte- 
grand tends to a constant corresponding to the £ = values of 
e(0) > 1, yu(0) ~ 1, and k(0) = 0. Thus in this d — ► oo limit 
E(d)/A oc d~ 3 and F(d)/A oc d~ 4 and the force is always 
repulsive, since essentially only e(0) matters. This analysis 
shows that it is crucial to employ the correct limiting values 
of e(i£), fi(i£,), as £ — > oo and £ — > 0, since these val- 

ues determine the behavior of E(d)/A in the limit d — ► and 
d — > oo respectively. 

In this work we have extended the Lifshitz theory to calcu- 
late the Casimir force by including chirality terms for the first 
time. We have shown that the chirality, if strong enough, is of 
critical importance in producing nanolevitations under realis- 
tic frequency-dependence and correct limiting values of e(u>) 
and fi(uj). Note, the previous calculations claiming repulsive 
Casimir force between metamaterials separated by vacuum 
have been achieved at the expense of nonrealistic frequency- 
dependence and/or limiting values of e(ui) and n(ui). Thus, 
chiral metamaterials might possibly be the main candidates to 
achieve experimentally the goal of Casimir repulsion, which 
might open up many opportunities for application. 
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